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Abstract 

The magnetic domain structure and magnetic properties of a ferromagnetic (Ga,Mn)As epilayer 
with perpendicular magnetic easy-axis are investigated. We show that, despite strong hystere- 
sis, domain theory at thermodynamical equilibrium can be used to determine the micromagnetic 
parameters. Combining magneto-optical Kerr microscopy magnetometry and ferromagnetic reso- 
nance measurements, we obtain the characteristic parameter for magnetic domains A c , the domain 
wall width and specific energy, and the spin stiffness constant as a function of temperature. The 
nucleation barrier for magnetization reversal and the Walker breakdown velocity for field-driven 
domain wall propagation are also estimated. 

PACS numbers: 75.60.-d, 75.50.Pp, 75.30.Gw 



1 



In thin layers based on ferromagnetic (FM) metallic or semiconducting materials do- 
main wall (DW) propagation is the subject of extensive research with the ultimate aim 
of information storage and/or transport by domain walls.- i2«MiIL£aL§, In particular, the FM 
semiconductor (Ga,Mn)As has received a lot of attention recently since the critical current 
for DW propagation was shown to be two orders of magnitude smaller than in metallic 
systems.- Both current- driven and field-driven DW propagation have been recently inves- 
tigated in (Ga,Mn)As layers and stripes.-^ In this context, it is important to determine 
experimentally the micromagnetic parameters such as the DW width and spin stiffness con- 
stant since they play a crucial role in DW dynamics. In the hydrodynamic regime for 
field-driven DW propagation, the DW mobility and the Walker breakdown velocity are pro- 
portional to the DW width,- which depends on the spin stiffness constant and the magnetic 
anisotropy constants.™ Spin stiffness also governs the exchange length, which characterizes 
the Bloch line width and sets an upper limit for the size of superparamagnetic particles. 
Moreover, the determination of the spin stiffness constant also gives access to the J p d ex- 
change parameter characterizing the p — d exchange interaction between localized Mn spins 
and itinerant carriers (holes) The experimental determination of the spin stiffness is 

not straightforward. It was extracted from spin wave resonances despite the discrepancy 
between theoretical and experimental dispersion curves.— 1 ^ The obtained value was larger 
than the predicted ones. In this paper, we determine the spin stiffness constant and the 
other micromagnetic parameters using combined experimental techniques and analysing the 
magnetic domain structure and the hysteresis cycles in the framework of domain theory. 

Owing to the complex structure of the valence band in zinc-blende (Ga,Mn)As, the di- 
rection of the magnetic easy-axis can be set parallel or perpendicular to the plane by tuning 
the layer strain and carrier concentration.— iI£iiliiiLi£ i n the case of strong perpendicular 
anisotropy, the ground state of the FM layer at thermodynamical equilibrium corresponds 
to the self-organized periodic stripe or bubble pattern of the up- and down-magnetized 
domains. 1 - Domain theory derived from the micromagnetic equations for FM films provides 
a quantitative description of the equilibrium domain patterns. Long-range self-organization 
results from the competition between the DW energy and the magnetostatic energy. The do- 
main pattern is controlled by a single parameter A c = i c jd = a/2K<id, where £ c is a character- 
istic length equal to the ratio of the specific wall energy a and the magnetostatic energy with 
a = 4V AK U and = fioM^/2. A is the spin stiffness constant, K u the uniaxial anisotropy 
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constant, M s the saturation magnetization, and d the sample thickness. Experimentally, the 
parameter A c can be extracted from the zero-field period of domain array. This method was 
used by T. Dietl et al. to determine A c for perpendicular-axis (Ga,Mn)As layers showing 
magnetic domains. 1 - However, the question arises whether the observed domain structures 
correspond to the equilibrium state. Actually, for perpendicular-axis (Ga,Mn)As the hys- 
teresis cycle is square,— 1 ^ thereby indicating that metastable states and/or DW pinning play 
an important role in the field dependence of the magnetization. Defect-assisted nucleation 
processes and DW pinning were recently investigated using magneto-optical imaging.— i22j2I 
In particular it was reported that the equilibrium demagnetized stripe or bubble state at 
zero applied field cannot be reached using AC demagnetization. 

In this paper we show that domain theory can be applied to determine upper and lower 
boundaries for the parameter A c even in the presence of metastability and DW pinning. 
The upper boundary is obtained by comparing the domain width close to magnetization 
saturation with the predicted equilibrium width. The lower boundary is obtained from the 
comparison of major and minor hysteresis loops with the equilibrium magnetization curve. 
Using the saturation magnetization and the anisotropy constants obtained by supercon- 
ducting quantum interference device (SQUID) magnetometry and ferromagnetic resonance 
(FMR) experiments, respectively, we obtain the temperature dependence of the lower and 
upper boundaries for the micromagnetic parameters: the specific DW energy, the DW width, 
and the spin stiffness constant, which are determined within a narrow range. 

The sample was prepared by molecular beam epitaxy. It consists of a 50 nm thick 
Gao.93Mno.07As layer grown on a Gao.902lrio.o9sAs relaxed buffer layer deposited on a semi- 
insulating (001) GaAs substrate. After post-growth annealing the Curie temperature Tq 
was 130 K. The magnetic domain structure was investigated using magneto-optical Kerr 
(MOKE) microscopy. The hysteresis cycle was obtained from the average intensity of MOKE 
images as a function of the applied field. Experimental details can be found elsewhere^. 

The anisotropy constants were obtained from FMR measurements. The FMR spectra 
were measured with X-band and Q-band spectrometers in the first derivative absorption 
mode with a modulation frequency of 100 kHz. To determine the anisotropy constants 
the spectra were measured for both the in-plane and out-of-plane orientation of the static 
magnetic field. From the resonance field positions we can deduce that the easy-axis of 
the sample is the along [001] direction for all temperatures T < Tc- For the in-plane 



3 



angular variation of the magnetic field an almost isotropic behavior is observed. Using 
the Smit-Beljers equation and the minimization of the free energy for different alignments 
of the magnetization^^ the magnetic anisotropy constants K 2 ±, K<$±, K 2 \\, and Km\ were 
obtained. They are plotted in Fig. ( dbj SIS db function of temperature. The small values 
of the in-plane anisotropy constants K 2 \\ and Km\ explain the almost isotropic in-plane 
behavior. In the following we will neglect them and consider the first term of the development 
of the anisotropy energy K u sin 2 9 with 9 the angle between the [001] direction and the 
magnetization. We obtain K u from the FMR results as K u = K 2 ± + K±±. The temperature 
dependence of K u is plotted in Fig. (a). The uniaxial anisotropy is characterized by the 
parameter Q = K u /Kd- From FMR and SQUID results we obtain a Q parameter ranging 
from 8.6 to 14 (Fig. (b)). 

Since Q > 1 domain theory for FM films with strong uniaxial anisotropy can be applied. 10 
In this framework the free energy of a stripe array F str is the sum of the DW energy, the 
Zeeman energy and intra- and inter-domain magnetic interaction terms. F str is a function 
of two variables: the reduced average magnetization along the direction of the applied field 
m = (M) /M s and the stripe period P. Minimization of F str provides the dependence of 
m and P on A c and on the reduced applied field h = H/M s . m(h) curves are shown in 
Fig. (a) for several values of A c . The stripe width of the minority phase (M opposite to 
H) decreases as h increases from zero up to the saturation field h s = hr m -+i) (Fig- (b))- It 
is worth noting that this width keeps a finite value W s when m — > 1 due to intra-domain 
magnetic interaction. The dependence of W s on A c is obtained by solving the equation 
27rA c = w 2 ln(l + w~ 2 ) + ln(l + w 2 ), where w s = W s /d. It is shown in Fig. (b). Let us 
note the quasi-exponential increase of the stripe width with A c , which means that domains 
cannot exist for A c ^> 1. 

In (Ga,Mn)As with perpendicular easy axis, DW pinning is strong enough with respect to 
inter-domain magnetic interaction to hinder the long-range ordering of magnetic domains. 17 
However, lamellar domains with a narrow distribution of their width are observed close 
to magnetization saturation, as shown in Fig. (c). This indicates that DW pinning is not 
sufficiently strong to prevent the formation of lamellar domains resulting from intra-domain 
magnetic interaction. An upper boundary for the parameter A c is obtained by comparing the 
prediction for the equilibrium domain width W s with the average width W of these lamellar 
domains. At T = 80 K, W is found equal to 1.7±0.3 /im, which corresponds to A c =1.3. 
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Since DW pinning impedes the reversible evolution of the domain width with the field, the 
width of lamellar domains should be larger than the predicted equilibrium stripe width. 
Moreover, owing to the limited spatial resolution (« 1/im), the measured width is probably 
larger than the actual one. Finally, since W s is an increasing function of A c , the value of A c 
deduced experimentally corresponds to an upper boundary for A c . The same analysis was 
repeated at different temperatures. The results are reported in Fig. (a). A small variation 
of \™ ax (1.24 to 1.7) is found as the temperature T varies from 12 to 120 K. 

A lower boundary for A c is determined by comparing the theoretical predictions for the 
magnetization curve m(h) with the hysteresis cycles obtained experimentally. Because of 
hysteresis, the m(H) curve predicted for thermodynamical equilibrium and without pinning 
should lie inside the experimental hysteresis cycle. A major hysteresis cycle at T=80 K is 
shown in Fig. (b) (black squares). The applied field is swept up to 60 mT, i.e., beyond the 
saturation field (~ 12 mT), which is much larger than the coercive field fioH c = 1.6 mT. 
When the field is swept down, the onset of magnetization reversal and therefore the width 
of the cycle is determined by defect-assisted nucleation (see below). In order to obtain a 
narrower cycle, we recorded a minor hysteresis loop (empty squares). It is obtained by 
sweeping the field up to fioH m = 6 mT, i.e., below the saturation field, in order to leave a 
few reversed domains in the sample. Hence, upon sweeping the field down, the nucleation 
stage is avoided. The shape of the hysteresis loop is determined by the dependence of the 
DW velocity on the applied field. Care is taken to use a low field ramp rate in order to 
allow for DW propagation at low velocity in the vicinity of the coercive field. The shape 
and width of the minor hysteresis loop does not depend on H m provided it is larger than the 
coercive field. The theoretical m(h) curves for stripes were calculated for a set of A c values 
(Fig. (a)). The saturation field increases for decreasing A c . Therefore, by selecting the m(h) 
curve which is tangent to the hysteresis cycle close to saturation, one can obtain a lower 
boundary for A c . This analysis can be further refined by considering that the saturation 
field calculated for stripes is lower than the theoretical collapse field, i.e., the field at which 
domains of the minority phase disappear. As m — > 1, lamellar domains keep a finite width 
but their length decreases until they reach the bubble circular shape and finally collapse.— 
The bubble collapse field h co u can be calculated from the free energy of an isolated bubble 
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with reverse magnetization^ 1 ^ 

F sb = 7rK d d 3 [2X c p +(h- l)p 2 + p 2 f(p)] , (1) 

where p = 2R/d is the reduced diameter of the bubble and f(p) = 1 + 
Ap{l - [(2k 2 - l)E(k) + (1 - k 2 )K(k)} £r 3 } /3tt with k 2 = p 2 /(l + p 2 ). E and K are 
complete elliptic integrals of the first and the second kind, respectively. F sb varies non- 
monotonously with p as shown in Fig. (a). Above the critical field h cr n the system re- 
mains is a metastable state until the collapse field is reached. h co u is the field at which 
dF sb /dp = d 2 F sb /dp 2 = 0. It increases with increasing A c . 

Figure (b) shows the the critical and collapse fields. The dashed curve shows the expected 
magnetization path from the stripe array curve to the bubble collapse field for the metastable 
domain pattern. The lower boundary for A c is taken as the value that makes this curve 
tangent to the minor hysteresis cycle close to saturation. At T = 80 K, A™ n =0.85 is 
obtained. This procedure was repeated at different temperatures. 

Finally we obtain the micromagnetic parameters. They are displayed in Figure . The 
specific wall energy a, the DW width nA = -n^J A/ K U) and the exchange spin stiffness 
constant A are derived from the A c values using the temperature dependence of M s and K u . 
a is found in the range 1.5-11.5 10~ 5 J m -2 , which is about 2 orders of magnitude smaller 
than for iron or cobalt films.— ^ This explains that propagating DWs easily skirt around 
pinning defects.— 1 ^ The DW width is found in the range 4-11 nm. This is larger than the 
mean distance between Mn ions (0.9 nm for 7% Mn) and smaller than the domain width, 
which validates the use of domain theory. The ratio of the DW parameter to the sample 
thickness A/d is larger than 0.025. Together with the large Q value, this ensures that the 
specific DW energy is very close to 4a/ AK u , as assumed in our analysis. 27 

From the micromagnetic parameters one can estimate the Walker breakdown velocity Vw 
for field-driven DW propagation in the flow regime.— Using V w = 7v /2 / u Ag( v /l + T/Q - 
1) ~ , ~ffioM s A/2 with 7 the gyromagnetic ratio, one finds Vw ~ 5.5 m s" 1 for the lower 
A-values and 15 to 9.6 m s _1 for the largest A-values for T between 12 and 80 K. These 
last Vw values are consistent with the largest velocities measured in the depinning regime 
for the same sample. 8 This indicates that the values of A c , a, A, and A are very likely the 
largest ones, i.e., those determined from the lamellar domain width. 

From the determination of A c , one also gains some insight into the magnetization reversal 
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process. We obtain the height of the intrinsic nucleation barrier. As seen in Fig. (a), below 
the critical field the energy of the system with a reverse bubble domain is lower than the 
energy of the saturated state. However the system can stay in the saturated state since 
there is an energy barrier AE for nucleation. Taking 0.85< A c <1.3 at T = 80 K yields 
280< AE/KT <1310 at h — h crit . The barrier height only weakly decreases with decreasing 
field. At the onset of magnetization reversal (h m —0.05) one still finds large values for AE, 
namely 194< AE/KT <660. This means that the system is highly metastable and that 
nucleation proceeds through local defects that strongly decrease the barrier height. This 
conclusion also holds at T = 120 K, i.e., close to Tc- 

The spin stiffness constant is smaller by one order of magnitude than the value determined 
by Goennenwein et al. from spin wave resonances in FMR spectra for (Ga,Mn)As with in- 
plane magnetization.— Let us note that a linear gradient of the magnetic properties had to 
be assumed in [ijjl in order to fit the results with theoretical spin wawe dispersion curves. 
Theoretical predictions for A were made in Ref. 13J for various carrier densities and for two 
values of the J p d exchange constant, taking a Mn concentration of 4.5%. For our sample 
we obtain an effective Mn concentration x e ff = 3.8% (we neglect the hole magnetization 
and assume M S (T — > 0) = x e ffN 5fiB, with iV the density of cations site). With this 
Mn concentration, our results for the spin stiffness are consistent with the extrapolated 

n 

theoretical estimations of [131] provided that a low value of the J p d exchange constant is used 
(Jpd < 50 meV nm 3 ). 

To summarize, we have applied domain theory to the case of a metastable system in 
order to obtain the ratio of DW specific energy and magnetic energy. Combining the study 
of domain structure, magnetization and anisotropy constants in a (Ga,Mn)As thin film 
with perpendicular magnetization, we have determined the micromagnetic parameters: the 
specific DW energy (1.5-11.5 10~ 5 J m -2 ), the DW width (4-11.3 nm) and the spin stiffness 
constant(0. 01-0.1 pJ m _1 ). From these results we estimate the nucleation barrier and the 
Walker velocity for DW propagation. 



7 



* e-mail: catherine.gourdon@insp.jussieu.fr; web site: nttp : / / www . insp . upm cTfr] 

1 A. Yamaguchi, T. Ono, S. Nasu, K. Miyake, K. Mibu, T. Shinjo, Phys. Rev. Lett. 92, 077205 
(2004). 

2 N. Vernier, D. A. Allwood, D. Atkinson, M. D. Cooke, and R. P. Cowburn, Europhys. Lett. 65, 
526 (2004). 

3 L. Thomas, M. Hayashi, X. Jiang, R. Moriya, C. Rettner, and S. P. Parkin, Nature 443, 197 
(2006) 

M. Hayashi, L. Thomas, Ya. B. Bazaliy, C. Rettner, R. Moriya, X. Jiang, and S. S. P. Parkin, 
Phys. Rev. Lett. 96, 197207 (2006). 

5 M. Yamanouchi, D. Chiba, F. Matsukura, and H. Ohno, Nature (London) 428, 539 (2004). 

6 M. Yamanouchi, D. Chiba, F. Matsukura, T. Dietl, and H. Ohno, Phys. Rev. Lett. 96, 096601 
(2006). 

7 H. X. Tang, R. K. Kawakami, D. D. Awschalom, and M. L. Roukes, Phys. Rev. B 74, 041310(R) 
(2006). 

8 A. Dourlat, V. Jeudy, L. Thevenard, A. Lemaitre, and C. Gourdon, J. Supercond. Nov. Magn. 
20, 453 (2007). 

9 A. P. Malozenoff, J. C. Slonczewski, Magnetic domain walls in bubble materials, Academic Press, 
New York, 1979. 

10 A. Hubert and R. Schafer, Magnetic domains, Springer, Berlin, 2000. 

11 T. Dietl, H. Ohno, and F. Matsukura, Phys. Rev. 63, 195205 (2001). 

12 T. Jungwirth, J. Sinova, J. Masek, J. Kucera, A. H. MacDonald, Rev. Mod. Phys. 78, 809 
(2006). 

13 J. Konig, T. Jungwirth, and A. H. MacDonald, Phys. Rev. B 64, 184423 (2001). 

14 S. T. B. Goennenwein, T. Graf, T. Wassner, M. S. Brandt, M. Stutzmann, J. B. Philipp, R. 
Gross, M. Krieger, K. Ziirn, and P. Ziemann, A. Koeder, S. Frank, W. Schoch, and A. Waag, 
Appl. Phys. Lett. 82,730 (2003). 

15 M. Sawicki, F. Matsukura, A. Idziaszek, T. Dietl, G. M. Schott, C. Ruester, C. Gould, G. 
Karczewski, G. Schmidt, and L. W. Molenkamp, Phys. Rev. B 70, 245325 (2004). 

16 L. Thevenard, L. Largeau, O. Mauguin, and A. Lemaitre, Appl. Phys. Lett. 87, 182506 (2005). 



8 



A. Dourlat, V. Jeudy, C. Testelin, F. Bernardot, K. Khazen, C. Gourdon, L. Thevenard, L. 
Largeau, O. Mauguin, A. Lemaitre, J. Appl. Phys. 102, 023913 (2007). 

L. Thevenard, L. Largeau, O. Mauguin, A. Lemaitre, K. Khazen, and H. J. von Bardeleben, 
Phys. Rev. B 75, 195218 (2007). 

T. Dietl, J. K6nig, A. H. MacDonald, Phys. Rev. B 64, 241201(R) (2001). 

L. Thevenard, L. Largeau, O. Mauguin, G. Patriarche, A Lemaitre, N. Vernier and J. Ferre, 

Phys. Rev. B 73, 195331 (2006). 

K. Y. Wang, A. W. Rushforth, V. A. Grant, R. P. Campion, K. W. Edmonds, C. R. Staddon, 
C. T. Foxon, B. L. Gallagher,J. Wunderlich, and D. A. Williams, J. Appl. Phys. 101,106101 
(2007). 

X. Liu, W. L. Lim, L. V. Titova, M. Dobrowolska, J. K. Furdyna, M. Kutrowski, and T. 

Wojtowicz, J. Appl. Phys. 98, 063904 (2005). 

J. Smit and H. G. Beljers, Philips Res. Rep. 10, 113 (1955). 

C. Gourdon, V. Jeudy, A. Cebers, Phys. Rev. Lett. 96, 087002 (2006). 

S. Huo, J. E. L. Bishop, and J. W. Tucker, J. Appl. Phys. 81, 5239 (1997). 

H. Chen and S. L. Whittenburg, J. Appl. Phys. 94, 5278 (2003). 

A. Hubert, J. Appl. Phys. 46, 2276 (1975). 



9 



FIGURE CAPTIONS 



FIG.l (a) Temperature dependence of the anisotropy contants K 2 ± (black squares), K 4 ± 
(black circles), K 2 \\ (open squares), K±\\ (open circles). K u = K 2 ± + K 4± (black diamonds, 
dashed line), (b) Temperature dependence of the magnetization (black squares, left scale) 
and of the Q parameter (black circles, right scale). 

FIG. 2 (a) Dependence of the reduced average magnetization m = (M) /M s on h = H/M s 
for a stripe array with A c = 0.6 (full curve), 0.8 (dashed), 1.2 (dash-dotted), (b) Stripe 
width as a function of the parameter A c for stripes of the minority phase near saturation 
(solid curve) and stripes at m — 0, H — (dashed curve), (c) Lamellar domain structure 
in (Ga,Mn)As near the saturation field at T=80 K. 

FIG. 3 Temperature dependence of the upper (black squares) and lower (empty squares) 
boundaries of : (a) the parameter A c , (b) the specific wall energy a, (c) the exchange spin 
stiffness A, and (d) the DW width nA. They are obtained from the lamellar domain width 
(upper values) and hysteresis cycle (lower values). 

FIG. 4 (a) Energy of an isolated bubble with reverse magnetization for several magnetic 
fields (A c =l). (b) Major (black squares) and minor (empty squares) hysteresis cycles ob- 
tained from Kerr microscopy images at T=80 K. The full curve is the m{h) curve for a stripe 
array with A c =0.85. The bubble critical and collapse fields for A c =0.85 are indicated as h crit 
and h co u, respectively. The dashed curve links the m(h) curve for stripes to the bubble 
collapse point (h co u, m = 1). It shows the expected magnetization path resulting from the 
existence of the collapse barrier. 
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FIGURE 2 
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FIGURE 3 
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